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Abstract 



The gaussian closure approximation previously used to study the growth 



0\ 

kinetics of the non-conserved 0(n) model is shown to be the zeroth-order ap- 
■ proximation in a well-defined sequence of approximations composing a more 



elaborate theory. This paper studies the effects of including the next non- 



C3 

a 

■ trivial correction in this sequence for the case n = 2. The scaling forms for 

a: 

the order-parameter and order-parameter squared correlation functions are 



determined for the physically interesting cases of the 0(2) model in two and 
three spatial dimensions. The post-gaussian versions of these quantities show 
improved agreement with simulations. Post-gaussian formulae for the defect 
density and the defect-defect correlation function g{x) are derived. As in the 
previous gaussian theory, the addition of fluctuations allows one to eliminate 
the unphysical divergence in g(x) at short scaled-distances. The non-trivial 
exponent A, governing the decay of order-parameter auto-correlations, is com- 
puted in this approximation both with and without fluctuations. 
PACS numbers: 05.70.Ln, 64.60.Cn, 64.60.My, 64.75. +g 
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I. INTRODUCTION 



In the gaussian closure approximation for phase-ordering kinetics the order-parameter 
field is expressed in terms of an auxiliary field which is assumed to obey gaussian statistics 
HUH- Recently, one of the authors extended this approach to treat more general statistics 
and applied the method successfully to the cases of conserved and non-conserved scalar 
fields HQ . This post-gaussian approach successfully eliminated several shortcomings of the 
gaussian theory. The present paper has two primary goals. The first goal is to generalize 
the post-gaussian theory to treat the non-conserved 0(n) symmetric model with continuous 
symmetry (n > 1). Systems with continuous symmetry have many physical realizations, 
including ferromagnets, superfluids, and liquid crystals j5]. We will focus on the case n — 2, 
where the defect-defect correlation function, g(x), has an unphysical divergence at short 
scaled-distances, x, in the gaussian theory. While the post-gaussian theory weakens this 
divergence, it is not eliminated. The second goal of this paper is then to show that the 
inclusion of fluctuations counter-balance and thereby eliminate the divergence in g(x) in the 
post-gaussian case. 

In a phase-ordering scenario, the dynamical evolution of the order parameter = 
ip(ri,ti) is not typically governed by a gaussian probability distribution, and analytical 
progress up to now has relied on relating t/>(1) to an auxiliary field assumed to be 

gaussian. The gaussian approximation has been very successful in treating the scaling inher- 
ent in the late-time behaviour of a host of growth kinetics problems 0. The gaussian theory 
predicts the now well-established result that for late times following a quench from the dis- 
ordered to the ordered phase the dynamics obey scaling and the system can be described in 
terms of a single growing length 

L(t) ~ t*, (1.1) 

which is characteristic of the spacing between defects at time t after the quench. <p is a growth 
exponent. In this regime the scaling form T for the equal-time {t\ = ti = t) order-parameter 
correlation function 
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CV(12) = (^(1) ■ tf»(2)> (1.2) 

can be accurately calculated within the theory. Here ?/>o is the magnitude ip = of the 
order-parameter in the ordered phase. The scaled length x is defined as x = r/L(t) with 
r = |r| = |r 2 — ri|. For the non-conserved models considered here the growth exponent 
is = 1/2, which is predicted by the theory and well established by experiments and 
simulations 0. The gaussian theory also makes quantitatively accurate predictions @,|||| 
for the exponent A governing the decay of the order-parameter auto-correlations, and defined 
by 

CV(0,M') ~ j^r-r for t»t'. (1.3) 

Finally, in addition to these accomplishments, the gaussian approximation is relatively easy 
to implement, and has straightforward generalizations to more complex systems. 

Despite these achievements there remain a few unresolved issues. The approximate 
nature of gaussian statistics was highlighted in the work of Blundell, Bray, and Sattler 



(BBS) (|T0| , pTl| where they computed, within the gaussian closure approximation, the two- 



point correlation function for the square of the order-parameter field: 

It is usual in comparing the theoretical scaling function T{x) with the data (both of which 
are relatively featureless) to rescale the length x to give the best fit. This rescaling reflects 
the uncertainty in the relation of the theoretical pre-factor of the power law growth ( |1 . 1| ) 
to the pre-factor determined in experiment and simulation [12]. By plotting l/(Cw,a + 1) 



against 1 — J 7 , hereafter referred to as a BBS plot, BBS were able to eliminate any adjustable 
fitting parameter and show that there were qualitative differences between the simulation 



results |] | and the predictions of existing gaussian theories (see Fig. y). These discrepancies 



indicated a need to go beyond the gaussian approximation. Another motivation for going 
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beyond the gaussian approximation is the ad hoc nature of the approximation itself and the 
desire to encompass it within a more general and systematic framework fT"3|| . 

A second problem with the gaussian closure approximation occurs when one examines 
defect correlations. In this paper we focus on point defects (n = d) whose density is defined 

as 

= $>« <*(ri - x^tO) (1-5) 

a 

where x a (ti) is the position at time t\ of the ath point defect, which has a topological charge 
q a . Defect-defect correlations 

0(12) = (p(l)p(2)> (1.6) 
at equal-times can be shown [O] to decompose into two parts 



g(r,t) = n (t)5(r)+g dd ( r ,t). (1.7) 

The quantity no(t) represents defect self-correlations and is just the total unsigned number 
density of defects at time t. We will be primarily concerned here with the defect-defect 
correlation function gdd{^,t) which measures the correlations between different defects. The 
conservation law 

Jd n rg(r,t)=0 (1.8) 

relates n (t) and gdd(*, t) through 

n (t) = - [ d n r g M {r,t). (1.9) 



In the scaling regime it can be shown |TJ| that n (t) ~ L n (t) and that ^(r,t) has the 
form 

9dd(r,t) = J2^j9( x ) (1-10) 

where g(x) is a universal scaling function. While the form for g(x) obtained from the 
gaussian theory for the two-dimensional 0(2) model is in good agreement with simulations 



lq| and experiments [jTIJ at large scaled-distances, there are qualitative differences in the 
short scaled-distance behaviour. The gaussian theory |TJ|] predicts a divergence in g(x) at 
small x while experiments and simulations have g(x) approaching zero at the origin (see Fig. 



This paper addresses these issues for the non-conserved 0(n) model by using a well- 
defined sequence of approximations for the probability distribution of the auxiliary field rh, 
which reduces to a gaussian distribution at lowest order. The theory presented here treats 
in detail the next non-trivial term in the sequence. This post-gaussian approach has been 
successfully applied to both the cases of a non-conserved and a conserved ||] scalar order 
parameter. While, as expected, the non-conserved scalar theory predicts little change in the 
form for J 7 , the BBS plot shows a marked improvement over the gaussian theory when com- 
pared with simulations. This encourages one extend the post-gaussian theory to the 0(n) 
case. The key result of this paper is that J 7 , C^a, no(t), g(x) and A can all be extracted using 
non-gaussian statistics. This is non-trivial since these quantities have universal forms in the 
scaling regime. This is connected to the fact that the probability distribution governing the 
auxiliary field in the gaussian case has a fixed-point form determined by the solution of an 
eigenvalue problem. In the post-gaussian case the determination of the fixed-point form for 
the probability distribution requires the solution of a double eigenvalue problem. This paper 
focuses on the 0(2) model in two and three spatial dimensions where experimental and sim- 
ulation results are readily available. For the 0(2) model in the post-gaussian approximation 
we find that the form for T changes little from the gaussian results and, as in the scalar case, 
the agreement of the BBS plot with simulations is improved. This improvement includes 
a dependence on dimensionality seen in the simulations, but not exhibited by the gaussian 
theory. There are some problems related to the negativity of C^i at small x but these are 
seen to be a consequence of the manner in which we defined the sequence of approximations. 
The exponent A is now in poorer agreement with the simulation data. The divergence in 
g(x) is found to be weaker in the post-gaussian theory than in the gaussian theory, but it is 
not completely eliminated. 



We recently showed that this divergence can be eliminated in the gaussian case [17| if 
one includes the effects of fluctuations about the ordering field. The origin of the divergence 
is the appearance of non-analytic terms in the small- x expansion of the auxiliary field cor- 
relation function /. One removes the divergence by eliminating the non-analytic terms in / 
through a careful treatment of the fluctuations. This development is theoretically pleasing 
since one expects the auxiliary field correlation function to be well-behaved. The question 
remains whether this post-gaussian scheme can be smoothly generalized to include these 
fluctuations. We answer this question in the affirmative and see that the post-gaussian 
theory with fluctuations is a rather natural generalization of the gaussian theory. Again, 
the divergence in g(x) is removed and the magnitude of g(0) is reduced, bringing the post- 



gaussian theory into better agreement with the simulation results |15j than the gaussian 
theory (see Fig. [|). 

The first part of this paper is mainly devoted to developing the results for the post- 
gaussian theory. Later, after we present the post-gaussian results for g(x), we will discuss 
the role of fluctuations in detail. Section II reviews the gaussian 0(n) model and the sepa- 
ration of the equation of motion into an equation for the ordering field and an equation for 
the fluctuations. The main results of this paper are contained in Section III which, after 
introducing the post-gaussian probability distribution and general formulae for calculating 
post-gaussian averages, derives post-gaussian expressions for J 7 , and the equations of 
motion. Defect correlations are discussed, leading into Section IV which considers the inclu- 
sion of fluctuations in the post-gaussian theory. The results of numerical analysis of the new, 
post-gaussian eigenvalue problem are presented in Section V. We conclude with a summary 
and discussion. 



II. MODEL 
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A. Preliminaries 



We consider the 0(n) model, which describes the dynamics of a non-conserved, n- 
component order-parameter field if>(l) = (0i(l), ■ • ■ , if) n (l)). To begin we will work with 
general n; however, later we will focus on the interesting case n — 2. As in previous work in 
this area ||, the dynamics are modeled using a time-dependent Ginzburg-Landau equation 

dt y ' 

We assume that the quench is to zero temperature where the usual noise term on the right- 
hand side is zero flIB| . T is a kinetic coefficient and F[if>] is the free-energy, assumed to be 
of the form 

F$] = J >r(||V$ a + V[if)}) (2.2) 

where the potential V[ip] is chosen to have 0(n) symmetry and a degenerate ground state 
with \if>\ = if) |I9| . With a suitable redefinition of the time and space scales the coefficients 
T and c can be set to one and fl2.1| ) can be written as 

It is believed that our final results are independent of the exact nature of the initial state, 
provided it is a disordered state. 

The evolution induced by (|2.3| ) causes if) to order and assume a distribution that is far 
from gaussian. It is by now standard to introduce a mapping between the physical field if) 
and an auxiliary field m with more tractable statistics. We can decompose if) exactly as 

■0 = a[m) + u. (2.4) 

The utility of this decomposition lies in our ability to create a consistent theory with the 
mapping a chosen to reflect the defect structure in the problem, and the fluctuation correc- 
tion u constructed to be small at late-times. The precise statistics satisfied by the fields m 
and u will be specified below. 



The mapping giving a as a function of fh is chosen to incorporate the dominant defect 
structure in the late-stage ordering kinetics. We assume that a satisfies the Euler-Lagrange 
equation for the free-energy given by Q2.2Q with the spatial coordinate replaced by the 
auxiliary field: 

V>1m] = (2.5) 



The defects are then the non-uniform solutions of (|2.5| ) which match on to the uniform 
solution at infinity Since we expect only the lowest-energy defects, having unit topological 
charge, will survive to late-times the relevant solutions to ( p.5|) will be of the form 



a[fh] = A(m)rh (2-6) 

where m = \m\ and m = fh/m. Thus the interpretation of rh is that its magnitude represents 
the distance away from a defect core and its orientation indicates the direction to the defect 
core. We expect m, away from the defect cores, to grow as L in the late-time scaling regime. 



Inserting ( |2.6|) into ( |2.5| ) gives an equation for A, 

T) — 1 

V 2 m A-?- ir A-V'[A)=0, (2.7) 

where the prime indicates a derivative with respect to A. The boundary conditions are 
A(0) = 0, A(oo) = i[)q. An analysis of ( PTT| ) for n > 1 and large m yields 



A(m) = ip 



K 

!- — + ■• 



(2.8) 



where k = (n — l)/V"[i/)o\ > 0. The algebraic relaxation of the order-parameter to its ordered 
value is a distinct feature of the 0(n) model for n > 1 which must be carefully treated in 
the evaluation of some averages. In the scalar case (n — 1) tp relaxes exponentially to ipo 
away from the defects. 

B. Separation of the equations of motion 



In [F7| it was shown how one can quite generally separate the equation of motion ( |2.3p 
into an equation for the ordering field a and an equation for the fluctuating field u. One 
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finds that a satisfies the equation of motion 

B = (2.9) 

with 

B = d t B- V 2 a + V 2 m a - 0. (2. 10) 

is chosen so that u is small in the scaling regime and thus represents a fluctuation. The 
equation of motion for u is 

on • 

-l = ^ Ur -W il [a]u 1 -Q t (2.11) 

where a sum over the index j is assumed and, to leading order in 1/L, 

W ij [a]=q 2 a i a j (2.12) 

which is purely longitudinal, q 2 . = V"[ip ] > 0. 

If we set G equal to zero in ( |2.10p we obtain the equation used previously to determine 



the a correlations 0. This choice decouples a and u. The equation for u would then 
separate into a (massless |2"0]] ) diffusion equation for the transverse piece ut and an equation 
for the longitudinal piece ul with a mass term —q 2 UL- However, as was seen in the gaussian 
theory pjli], the equation for a would necessarily lead to non-analytic behaviour in the 



short scaled-distance expansion for the normalized auxiliary field correlation function, f(x), 
and would ultimately lead to an unphysical divergence in g(x) at small x. We must choose 
so that our auxiliary field (or fields, in the post-gaussian theory) is analytic for small x. 
The possible forms we can use for 9 are discussed in detail in [0. The key ideas are that 
(i) is chosen to be a function of rh only. This choice means that a satisfies a closed 
equation, while u is slaved by rh. This highlights the fact that, since we are working at zero 



temperature, it is and not temperature which is driving the fluctuations |2T 
(ii) must be odd under m — > —rh. 
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(iii) must scale as 0(L~ 2 ) in the scaling regime if it is to compensate for the terms in 
the equation of motion which lead to the non-analyticities in /. This will also allow us to 
treat fluctuation since it will imply u ~ L . 

It is sufficient for our purposes to consider 

e = ^, (2.i3) 

where Uq is a constant to be determined. This form will allow us to construct f(x) to be 
analytic through terms of 0(x 2 ). It was shown in how (|2.13|) is the first term is a series 



for which can be used to enforce analyticity at 0(x ) and beyond. 



C. The gaussian approximation 

To complete the definition of the model one must specify the form of the probability 
distribution for the auxiliary field m. Forcing a to satisfy the exact equation of motion ( |2.9| ) 
is tantamount to solving the problem exactly, and will determine a probability distribution 
for m which is complicated and extremely difficult for purposes of computation. Progress 
can be made if one imposes the weaker constraint 

(B(l) ■ <7(2)>o = 0. (2.14) 

This equation allows one to insure that B(l) is reasonably small at late-times but gives one 
the flexibility to choose a suitable probability distribution. The simplest choice is a gaussian 
probability distribution for m: 

P[m) = P [m) = Ne~ Ko[ ™ ] (2.15) 

with M a normalization constant and 

K [m] = ^J d\d1 C \l2)m{l) ■ m(2). (2.16) 

Co(12), the correlation function for the auxiliary field, is explicitly defined through 
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(m i (l)m i (2)) = Sa C (12). 



(2.17) 



Here we have used (• • -} to indicate an average using the gaussian distribution Po[^]- Later, 
when we consider post-gaussian statistics, we will use (• • •) to denote an average using the 
post-gaussian distribution P[m]. The system is assumed to be statistically isotropic and 
homogeneous so Cq(12) is invariant under interchange of its spatial indices. For future 
reference we also define the one-point correlation function 



So(l) = Co(ll) 



(2.18) 



and the normalized auxiliary field correlation function 

0,(12) 



/(12) 



So(12) 



(2.19) 



with 5 (12) = y<So(l)<So(2). As discussed above it is expected that both C and So grow as 
L 2 at late times. 

In the gaussian theory the relationship between the auxiliary function / and the observ- 
able functions C^, C^p. and g can be derived without reference to the dynamics contained 
in ( |2.14j ). Using Q2.4|), ( |2.6[ ) and ( |2.8j ) can be written to leading order in 1/L as 



CV(12)=^(7n(l).m(2)) 



(2.20) 



with 



2vr 



1 n+ 1 
2' 2 



1 1 n + 2 



2'2' 



;/ 2 



(2.21^ 



where S is the beta function and F is the hypergeometric function |2,22]. For later conve- 
nience, we define the short-forms 



F 1 = F 



2'2' 2 



(2.22) 



and denote 
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pi _ dF\ „ _ d 2 F x 
1 " df ' 1 ~ df 2 ' ' ' ' 

The quantity C^a ( |1.4| ), as obtained by BBS [Kj for n > 1, is given by 



^2.23) 



1 1- -■ f 2 

2 . 



1. 



(2.24) 



For later notational simplicity we shall write 



1 1--- f 

1,1, 2 ,7 



(2.25) 



and, as with i^, use the prime notation for differentiation with respect to /. As pointed 
out in ]TI , one can directly test the assumption of gaussian statistics, independent of the 
spatial form of /, by plotting 1/(C^2 + 1) against \ — T. This is done for the gaussian 0(2) 
model in Fig. |3| and the discrepancy between the gaussian theory and the simulation data 
is evident. Since this discrepancy is due to the choice of gaussian statistics and not due to 
the method used to determine / it strongly suggests that an improved choice of probability 
distribution is needed. 

Within the gaussian theory, g(x) is given by fil l 



with 



h = _l_ d J_ 

2tc dx 



(2.26) 



(2.27) 



and 7 = 1/ a/1 — p. The defect density is 



n {t) 



with 



2"W 2 r(l + n/2) 



= -([Vmf)„. 



D 



(2) 



nS (t) 



re/2 



(2.28) 



n 



(2.29) 



In the gaussian theory without fluctuations = 1. The expressions ( [2.26 ) and Q2.28 ) 
satisfy the conservation law ( |1.9|) . 
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The constraint ( 2.14T) determines the time evolution of the two-point order-parameter 



correlations. We use ( 2.14| ) to determine / and J 7 , which are related through ( |2.21| ). Knowl- 
edge of / allows us then to determine C^p. and g. To simplify the discussion we restrict 



ourselves initially to the case = 0, considered in ||. The case when has the form ( |2.13| ) 



was addressed in |T7|], and later we will discuss the inclusion of fluctuations in some detail. 
After some manipulation ( |2.14|) becomes 

d tl F(12) - V^(12) - ^—fd f F(12) = 0. (2.30) 

Here we have used the short-hand notation dfj 7 = dJ-'/df. For t\ ^> t 2 both JF and / are 
small. In this limit (|2.30|) becomes a linear equation for T and, with the definition 

L\t) = = At (2.31) 

2\x 

for the scaling length L f23|, A can be determined as |2|||: 



\ = d-— . (2.32) 

A/jL 

To examine the equal-time order-parameter correlations in the late-time scaling regime we 
set ti = t 2 = t and write (|2.30| ) in terms of the scaled distance x. To leading order in 1/L 



we have 



This is a non-linear eigenvalue problem with eigenvalue /i, entering via the definition (|2.31 



for the scaling length, fi is determined by numerically matching the analytically determined 
behaviour at small-x onto the analytically determined behaviour at large-x. For large x 
both T and / are small and ( |2.33| ) can be linearized. In this regime the physical solution to 
O) is 

T ~ x d - 2X e~ x2/2 . (2.34) 
These results are valid for arbitrary n > 1. Our focus, however, is on the 0(2) model where 



there are known qualitative discrepancies with simulation and experimental data fl4|--|16 
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With this in mind, we now examine the small- x behaviour of the scaling equation ( |2.33| ) for 
the case n — 2. For small- x (|2.33|) admits the following general expansion for /: 



f = l + f2X 2 



l + ^-[.l. + C) 

mx 



In a; 



1 + ^(1 O 

mx 



In x 



0(x 6 ). (2.35) 



Non-analyticities in / appear as a result of the non-zero K 2 and K4 coefficients multiplying 
factors of 1/lnx. 

The coefficients of the expansion Q2.35| ) can be determined by examining Q2.33| ) order- 
by-order at small x. Balancing terms at 0(lnx) gives 



h 



71 



[2.36) 



At (9(1) one has 



d' 



^2.37) 



The non-zero K% coefficient is particularly important since it is responsible for the divergence 



of the defect-defect correlation function at small x, as one can see by examining (|2.26|) at 
small-x where one has, for n = d = 2, 



f 2 K 2 



2tt 2 x 2 (lnx) 



+ ■ 



(2.38) 



III. THE POST-GAUSSIAN THEORY 

A. Beyond the gaussian approximation 

Any scheme that claims to improve upon the gaussian approximation should include the 
following elements: it should be systematic, with the gaussian approximation being 'lowest 
order', it should be fairly easy computationally, at least in the early stages, and it should be 
converging to the results of experiment and simulation the farther one goes in the scheme. 
We now present an approach which generalizes the previous gaussian theory and satisfies 
the first two criteria completely, and the third criterion in all areas save one. 
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Since any functional can be written as a sum over generalized Hermite functional poly- 
nomials we write the probability distribution for the auxiliary field in the form || 

oo 

P[m] = P [m]J2 E aj{i 1 ,---,ij;l---J)Hj{i 1 ,---,i J ;l---J) (3.1) 

J=0 

where the indices i\ ■ ■ ■ i j each range from 1 to n. Integration over repeated spatial and 
temporal variables is assumed. Po[^] is the gaussian distribution. The Hj are the Hermite 
functional polynomials 

Hj{i x , ■ ■ ■ , ij; 1 • • • J) = (-l)VToH _ e -*>[»] (3.2) 

dmi^l) ■ ■ ■dmi J (J) 

which form a complete orthogonal set, spanning a space containing the 0(n) symmetric 
functionals [[J. The functions aj, along with Cq, are determined by the symmetry of the 
problem and the series of constraints 

(B(l) ■ *(2)> = (3.3) 
(5,(1)^(2)^(3)^(4)) = (3.4) 



In the gaussian theory (J = 0) only the first constraint (|3.3|) was necessary to completely 



determine the dynamics of Co- As one does computations at higher J more of these con- 
straints are necessary to determine all the a,j. At each J one has a systematic approach to 
the problem and one can, in principle, calculate to any order in J. One expects that by 
enforcing more constraints on B one will be satisfying the exact equation of motion ( |2.9|) 
more stringently. At the same time one is developing a more accurate expression for P{m\ 
as more Hermite polynomials are included. In this sense then the theory is expected to 
improve as one calculates to higher J. Finally, the use of Hermite polynomials allows one 
to straightforwardly express post-gaussian averages in terms of easily computable gaussian 
averages. 

We will work to order J = 2, which we call the first post-gaussian approximation, so we 
will make the choices 
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a = 1 (3.5) 

ai(i;l) = (3.6) 

a 2 (iJ;12) = 6 ij A 2 {12) (3.7) 

Oj(ii, • • • , ij; 1 • • • J) = for J > 2. (3.8) 



Here A 2 (12) is a scalar function, symmetric in its arguments. The condition on ao insures 
that the theory reduces to the correct gaussian limit at lowest order and normalizes the 
probability distribution, provided Pq [m] is normalized. The condition on a\ reflects the fact 
that there are no external fields and P[fh] = P[—m\. The choice for a 2 follows from the 
0{n) symmetry, and considerations of isotropy and homogeneity. 

B. General results for post-gaussian averages 

In order to calculate physical quantities like Cf in the post-gaussian approximation we 
must be able to express one- and two-point averages like (0(1)) and (0(l)x(2)), where 
4> and x are functions of m, in terms of related gaussian averages. The calculation for the 
case with no spatial gradients in the average is presented below. The important case where 
there are spatial gradients in the average is slightly more involved, but the results for the 
0(n) model are straightforward generalizations of the results for the scalar case, presented 
in ||. The two-point average (0(l)x(2)) can be related to gaussian averages by using the 
definitions (|3.1|) and ( |3.2|) and doing a few integrations by parts. The result is 

(0(i)x(2)> = t E • • • , is, i • • • J) L ^t (1) A (2 \n ) • ( 3 - 9 ) 

J^o h%j \5m h (l) ■ ■ ■ 5 mij (J) / Q 

The barred quantities are integrated over in this expression, but they will only contribute if 
they take the values 1 or 2. In the first post-gaussian approximation ( |3.5| - |3.8| ) the two-point 
average ( |3.9|) is 

(0(1)X(2)) = (0(l)x(2)) o + A 2 (ll)(V^(l)x(2)) 

+A 2 (22)(0(l)V^ X (2))o + 2A 2 (12)(V m 0(l) ■ V mX (2)) , (3.10) 
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with no integrations over the spatio-temporal variables. If we now proceed under the as- 
sumption that A 2 (ll) and A 2 (22) are non-zero we will only be forced later to set them to 
zero so that the long-distance behaviour of the theory is physical. For notational convenience 
we set these quantities to zero at the outset. Equation ( |3.10| ) can be further simplified by 
the introducing the operator 

<5 < 12) = 1 + 2 - 42 < 12) ^ (3U) 
since, for a general gaussian average (0(l)x(2))o, the following identity holds: 

(V m 0(l) • V mX (2)) = (0(l)x(2))o. (3.12) 

The two-point post-gaussian average then can be compactly written 

(0(1) X (2)) = G(12) (0(l)x(2))o. (3.13) 

The operator notation illuminates the close relation between gaussian and post-gaussian 
averages in the first post-gaussian approximation. By setting %(2) = 1 in ( |3.10| ) one obtains 
the formula 

(0(1)) = (0(l)) o (3.14) 
for the one-point average in the first post-gaussian approximation. 



C. Cw, and C^2 in the first post-gaussian approximation 



The post-gaussian analogue of ( |2.20| ) for the order-parameter correlation function can 
be straightforwardly calculated using (|3.13|) . We write 



^(12)=V 2 G(12)(m(l)-m(2)> 
= ^(12) 



(3.15) 



where now 
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n 



T = —B 

2tt 



1 n + 1 

2' 2 



[/Fl + 2 ^(F 1 + /F 1 ')]. 



(3.16) 



The post-gaussian quantity g is denned as 



0(12) = 



A 2 (12) 



(3.17) 



in analogy with the definition for / (|2.19|) . Although / retains many features from the 
gaussian theory and g carries much information about the post-gaussian corrections, later 
we will find that / and g influence each other strongly, and must be determined together 
using the appropriate constraint equations. 



Equation ( |3.13|) can also be used to determine the post-gaussian G^p. from the gaussian 



C^2. We have 



Gp2 = F 2 + 2gF' 2 - 1 for n > 2. 



(3.18) 



Unlike the situation in the gaussian theory, we now need to know the specific forms for 
/ and g in order to create a BBS plot. Knowledge of the adjustable parameter in the length 
scale is, however, still unnecessary. Since the forms for / and g are needed in the BBS plot, 
one hopes that the shape of the BBS curve will now depend on spatial dimensionality. We 
should note at this point that if we had maintained A 2 (ll) ^ we would now be forced to 
set it equal to zero so that both and C^p, decay to zero at large distances, as is expected 
physically. 



D. Equations of motion 

To determine the unknown functions / and g we will use the constraints (|3.3j ) and ( p.4|) , 
enforced to keep B small. For now, we will work with = to make the post-gaussian 
analysis more transparent. 

From our experience with the scalar case 0,0, we know that the constraint equation 
( ^.3[ ) evaluated at a single space-time point, 

(B(l) ■ a(l)) = 0, (3.19) 
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contains information about the short-distance behaviour of the theory that allows us later 
to simplify the constraint equations for / and g. Written in full, ( |3.19| ) becomes 

(3.20) 



\d tl (a\l)) ~ <*(1) • V?*(l)> + {ail) ■ V^x(l)> = 0. 



Evaluating fl3.20| ) to leading order, we have 



d. 



(2) 



1 for n > 2 



(3.21) 



where 



4 2) = 4 2) + 24 2) 

4 2) = lim-V^ (12) 
4 2) = lim-V?A 2 (12). 



(3.22) 
(3.23) 
(3.24) 



The definition ( |3.23j ) for is equivalent to (|2.29|) if one considers ( p.29|) to be an average 
using a gaussian probability distribution with the post-gaussian C*o(12) as its correlation 

(2) 

function. The new constant A 2 will turn out to be determined as part of an eigenvalue 
problem. Since Sq is a manifestly positive quantity, ( |3 . 2 1| ) also provides an important 
upper bound on the eigenvalue A 2 . We must have A 2 < 1/2 for n > 2. 
We now examine the constraint ( |3.3| ) written out in full, 

d t MV ■ *(2)> " V?(a(l) • a(2)) + ([V^a(l)] ■ 9(2)) = 0, (3.25) 

which in the gaussian case is the sole equation needed to determine the dynamics of C$. 
The evaluation of ( |3.25| ) involves a straightforward application of (|3.13|) to the appropriate 
gaussian averages. At late times the leading order result is 



9^(12) - V?^(12) - 



n 



2tt5 (1) 



B 2 



1 n + 1 
2' 2 



[(/ + 2g)(F 1 + fF[) + 2^(2F{ + fF[')} = 0. 

(3.26) 



For t\ ^> t 2 J 7 , f and g are small and ( |3.26| ) again becomes a linear equation for T . The 
relation between A and /i is ( |2.32p , unchanged from the gaussian theory, and the definition 
of the scaling length L ( |2.31|) is retained. 
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As in the gaussian case we can write ( |3.26| ) as an equal-time scaling equation 

1 n + 1" 



4/i 



2' 2 



[(/ + 2g){F x + fF{) + 2gf(2F{ + fF{')} = 0. (3.27) 



Although the functions J 7 , f and g are related through ( |3.16|) and ( p. 27] ), the additional 
constraint equation (|3.4|) is needed to complete the theory and determine these functions 
separately. In (|3.4| ) we have a choice of how to contract the indices ijkl and 1234. It is 
important to note that the function g entering the first post-gaussian theory is a two-point 
quantity. Thus, unlike the usual case in perturbation theory, the first-order corrections to the 
gaussian theory will not require us to treat the difficult intricacies of four-point correlation 



functions. Therefore, in order to determine g we only need enforce (|3.4j) contracted to a 
two-point function. A non-trivial constraint is obtained by contracting the indices ijkl in 
pairs. There then remain two possible constraint equations: 

(B(l) ■ a{l) a{2) ■ a{2)} = (3.28) 
(B(l) ■ a{2) a{l) ■ a{2)} = 0. (3.29) 

Unless there exists some degeneracy one cannot, in the first post-gaussian approximation, 
satisfy relations ( |3.28|) and ( |3.29| ) simultaneously. Analysis of Q3.29| ) shows that for n = 2 



it produces an equation in which the eigenvalue A 2 does not appear |23] . This is contrary 



to our expectation, based on previous work 0,f|, that we have to solve a double eigenvalue 
problem in the post-gaussian theory. It is therefore clear that we should satisfy fl3.28p . 
Since we are enforcing ( |3.19| ) we may rewrite ( |3.28| ) as 

(.B(l) • <?(1)A(2)) = (3.30) 

where 

A(2)=a 2 (2)-^. (3.31) 

This computationally convenient form allows the calculation to proceed in a way similar to 
that in ||. Written in full, (|3~30|) becomes 
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([d h a(l)) ■ *(1)A(2)> - ([V?a(l)] ■ <?(1)A(2)> + ([V^(l)] • <?(1)A(2)> = 0. (3.32) 



To evaluate the post-gaussian averages in (|3.32|) one must use the formulae given in Section 



III.B and generalize the results in for post-gaussian averages containing spatial gradients 
of a scalar field m to the case of a vector field m. Two new gaussian averages specific to 
the 0(n) model must be calculated. After some algebra and rearrangement, ( |3.32|) reduces 
at late-times to an equal-time scaling equation for / and g. For n = 2 we have, at leading 
order in l/L, 

47r4 2) <7/ + 4fiVJV x g + ^{V x ff{l + Sgfr 2 ) = 0. (3.33) 

For n > 2 we have 

4nA^gf 2 F^ 2 + 4fiV x fV x g [2f 2 F 2 - (n - 2){F 2 - 1)] + 

/i(V,/) 2 [2f 2 F 2 - (n - 2)(F 2 - 1) + 2g(A 1 2 fF 2 + 2f 2 F' 2 - 

(n - 2)(Fi + 2( 7 2 - 2){F 2 - 1)//))] = 0. (3.34) 

Note that for n = 2, F 2 = 7 2 and ( p34l) reduces to ( p33|) . 

These equations, together with ( |3.16| ) and ( |3.27| ) form a complete set of relations that 



will be used to determine the functions f and g. There are two unspecified constants in 
this set of equations - /i and A 2 . We thus have a non-linear eigenvalue problem in which /z, 

(2) 

the eigenvalue familiar from the gaussian theory, and the new eigenvalue A 2 are determined 
by connecting the small- and large- a; behaviour of ( p.27|) and (|3.34|) . 



For large x (|3.27| ) reduces to a linear equation, as in the gaussian case, and once again 
leads to (|2.34|) . The form for the exponent d — 2A in (|2.34|) appears to be robust. An 



examination of equation (|3.34j) at large-x yields 

/ = -4<7- (3.35) 

Up to now these results have been valid for n > 1. Now, in examining the short scaled- 
distance properties, we will focus on the n = 2 case where there are logarithmic corrections. 
The generalization of fl2.35|) is 
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f(x) = 1 + f 2 X< 



mx 



O 



lnx 



1 + ^(1 O 



lnx 



+ 0(x 6 ) (3.36) 



g(x) = g 2 x 



l + r^-( 1 +0 



lnx 



+ 9aX 



1 + 7^ (10 

mx 



lnx 



0(x 6 ). (3.37) 



An examination of Q3.33D yields, at 0(x 2 ), 



7cA ( i ) g 2 + l x(f 2 ) 2 = 0. 



At C(x 2 /lnx) one has 



L 2 = 2K 2 . 



Examining ( p. 271 ) yields, at C(ln 



x 



h + 2^2 



7T 



while at 0(1), using ( |3.39| ) one has 

K,= 1 



7T 



2-d g 2 
2d h 



(3.38) 



(3.39) 



(3.40) 



(3.41) 



2 4/zd(/ 2 + %2) 

Equations Q3.38Q and ( |3.40|) determine / 2 and g 2 separately in terms of the eigenvalues /1 



and A^ 2) 



h 



71 



4/id 



A^d-y/A^>d(2 + A^d) 



I (2) 



( 2 ), 



(3.42) 



7T 



#2 



1 + A ( 2 2) d - \fA ( 2 >d(2 + A [ 2 'd) 



(2), 



1 (2) 



(3.43) 



8 fid 

We have assumed that A^ > and taken the negative square root in (ft.42 ) in order to 
render f 2 negative, as we expect / < 1 physically. Equations ( |3.39|) and (|3.41 ) determine 
the corrections to the leading-order behaviour. The small-x expansion of the scaling form 
( PIED for C r is 



C,^ (x) 



2^2 - h 1 , K 2 1 



2(f 2 ) 2 x 2 2f 2 x 2 \nx 







(x lnx) 



(3.44) 
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Note that the condition for to be positive near the origin, the behaviour expected 
physically, is 2g 2 — f 2 > which is true for 

4 2) > l/4d. (3.45) 

(2) 

We will find, however, that A 2 does not satisfy this lower bound, and the resulting nega- 
tivity of C^2 has to be interpreted carefully. 



E. Defect correlations 



The starting point for evaluating defect-defect correlations is to note [14| that the density 



for point defects ( |1.5|) can be rewritten in terms of the order-parameter field, 

p(l)=5[^(l)]det[V 1 ^(l)], (3.46) 



since the zeros of ip are the locations of the defects. From (|2.6 ) and (|2.7|) we see that ip ~ m 



near the defect cores so Q3.46j ) can be expressed in terms of the auxiliary field 



p(l) = 5[m{l)} det[Vim(l)]. (3.47) 

This form is convenient for evaluating the equal-time defect-defect correlations (|1.6| ) which 
separate as indicated in ( |l~7|) into a piece representing the defect self-correlations, 



n (t) = (5[m(l)]|det[V 1 m(l)]|), (3.48) 

and a piece representing correlations between different defects 

g dd (r,t) = (5[m(l)]5[m{2)} det[V x m(l)] det[V 2 m(2)]). (3.49) 

We sketch the calculation of gdd( r ,t) for the post-gaussian theory in the Appendix. We 
recover the scaling relation ( |1.10| ) with the post-gaussian scaling form 

x (h\ n ~ x (dh dh\ . 1X (h\ n ~ 2 hdh 

g{x)=n\[-\ + ) + (n ~ l)n\ - -— (3.50) 

\X I \OX OX I \X I X ox 



23 



where h has the same definition it had in the gaussian theory (|2.27 ) and h, which contains 
the new post-gaussian terms, is given by 



h 



7T \ dx 



dx 



(3.51) 



The defect density n (t) can be calculated directly by evaluating the post-gaussian average 
in ( 3.48Q . The derivation is similar to that given for the gaussian theory in |14| and, as in the 
gaussian theory, the absolute value of the determinant appearing in (|3.48| ) is a complication 
that has to be carefully treated. Alternatively, we can use the conservation law ( |1.9| ) to 
compute n (t), using the post-gaussian formulae ( |3.50| ) for g(x). Both approaches lead to 
the same result: 



n (t) 



2"W 2 r(l + n/2) 



o(2) 
°0 



n/2 



1 + n 



4 2) 



(2) 



(3.52) 



nS (t) 

Using the small- x expansions ( |3.36| ) and (|3.37|) for / and g we calculate the post-gaussian 
modification to the behaviour of g at small x for n = d = 2: 

{h + Ag 2 )K 2 1 



g( x ) 



+ ■ 



(3.53) 



2tt 2 x 2 (lnx) 2 

Thus, as in the gaussian theory, leading order non-analyticities in the small-x behaviour of 
the post-gaussian / and g are responsible for an unphysical divergence in g(x) at small- 
x. We will see that in the post-gaussian theory the divergence is weaker than it was in 
the gaussian theory; however, through the inclusion of fluctuations we can eliminate the 
divergence altogether. 



IV. THE POST-GAUSSIAN THEORY INCLUDING FLUCTUATIONS 

A. Analysis of the o degrees of freedom 

We will now show how the inclusion of fluctuations about the ordering field a can be 
used to eliminate the leading order non-analyticities in the small-x behaviour of / and g, 
thereby rendering g(0) finite. This has been done successfully for the gaussian theory in 
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|T7| . Fluctuations influence the equation of motion for a ( |2.9| ) through the non-zero G field 
( [2.131 ). The one-point equation ( |3.20|) , used to determine is modified by fluctuations to 



1 



~d t > 2 (l)) - ■ V?*(l)> + <*(1) ■ V^(l)) - 



^0 / _L> 



<a 2 (l)> = 0. 



(4.1) 



The only difference between ( |4.1| ) and ( |3.20| ) is the last term, which is 0(L 2 ). For n = 2 
the leading order contribution to ( |4.1| ) is 0(ln L/L 2 ) so the fluctuation term does not modify 

the result d = 1 for n = 2 in the post-gaussian theory. A similar thing happens in the 

(2) 

gaussian theory where the inclusion of fluctuations does not modify the relation Sq — 1 for 
n = 2 0. 

Fluctuations do, however, modify the two-point equation of motion (|3.25| ) determining 
the order-parameter correlations and lead to a new formula for A: 



A = d-f 

4/x 2 



(4.2) 



At equal-times in the scaling regime ( |3.25| ), including fluctuations, becomes 



n 



x ■ V X F + V l x T + u) Q T + — B 2 



1 n + 1 

2' 2 



[(/ + 2< ? )(F 1 + /F x ') + 2^/(2^ + /Ff)] = 0. 

(4.3) 



The final equation of constraint ( p.32|) is modified to 



([9 fl a(l)].a(l)A(2)>-([VMl)]'^(l)A(2)) 



LO 



L 2 (h 



■(a 2 (l)A(2)) + ([V^(l)]-a(l)A(2)) = 0, 

(4.4) 



however a calculation for n = 2 shows that at leading order in 1/L, this equation is un- 



changed from (|3.33j), even in the presence of fluctuations. Again, one has 



AirAfgf + 4/iVa/V^ + /i(V x ./) 2 (l + %gh 2 ) = 0. 



(4.5) 



The large- a; behaviour of the theory with fluctuations is essentially unchanged from the 
original post-gaussian theory. Equations ( |2.34| ) and ( |3.35p still hold except now A is given 
by ( |4.2|) . The expansions (|3.36|) and (|3.37|) are used to examine the small-x behaviour of 
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equations ( |4.3|) and ( |4.5| ). Once again f 2 and g 2 are given by ( |3.42|) and ( |3.43|) . The relation 
( p.39| ) still holds and from ( |4~3| ) one has 



loq = g 2 [2d + 



71 



71 



+ 2dK 2 (f 2 + 4g 2 ). 



(4.6) 



2fif 2 J 2fid 

At this point we insist that / and g are analytic for small-x at 0(x 2 ) (i.e. K 2 = L 2 = 0). 
This then fixes ujq in terms of /i and A 2 . 



u = g 2 \2d 



7T 



7T 



(4.7) 



2fif 2 J 2/jd 

The next order terms are 0(x A ). From (|4.3|) and ( |4.5|) we obtain two equations determining 
/ 4 and g A : 



(tt4 2) + ±nf 2 )g A + 4//(/ 2 - y 2 )/ 4 = 2^g 2 {f 2 f - tt^^ 2 / 2 



(2). 



(4.8) 



h + 2 ^ 4 = ^(/ 2 ) 2 + 3/ 2 ^ 2 + (/ ; ± 292 } 

J y 4 W 7 J y 4(d + 2) 



7T 



2-^o 

4/i 



(4.9) 



Similarly, at the next order two relations can be derived for K4 and L 4 , which are more 
involved. 

By using ujq to eliminate the terms at 0(x 2 / hxx) in / and g we have also managed to 
remove the divergence in g(x) at the origin. For small- 2 we now have 



7T Z V 2 



- - 2f 2 g 2 ) +o(±- 



(4.10) 



B. Analysis of fluctuation correlations 

Correlations in the fluctuation field u are completely determined within the theory. There 
are two types of equal-time averages that are of interest to us. The first describes cross- 
correlations between the a and u fields and is defined as 

C u0 {12) = (tf(n,t).*(r 2 ,t)). (4.11) 

The second describes correlations of the fluctuation field with itself and is given by 
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5 i:j C uu {12) = (uiir^u^t)). (4.12) 



As we will see later these quantities are closely related in the scaling regime. In [|I7| it was 
shown how one can form equations of motion for both C u0 and C uu by using the equations 
of motion ( f2~9[) and (|2.11|) for a and u. We can then use these equations to determine C u q 



and C uu explicitly if we make the additional assumption that u is a gaussian field coupled to 
the post-gaussian field fh. To effect this change we replace the gaussian functional e~ K °^ 
in equations ( |3.1|) and ( |3.2|) by the gaussian functional e~ Ko ^ m ' u ^ which is quadratic in both 



fh and u. Post-gaussian averages over fh are evaluated as before, while repeated use of the 
following identity 

( Ui (l)A[m,u\) = Jd3 C um (l3)(— j^A[m,u}) + J d3 C uu (l3)(j^A[rh,u]), (4.13) 
with 

<%C um (12) = ( Ui (l) mj (2)), (4.14) 

allows averages over u to be expressed in terms C uu , C um and post-gaussian averages over fh. 
We then evaluate the averages over fh, which can be expressed in terms of / and g obtained 
previously. Finally, we examine the equations of motion for the fluctuations in the late-time 
scaling regime and extract the scaling functions. The analysis closely follows that given in 
[T7| so here we report only the final results for the post-gaussian theory. 

In [^] it was shown that, as a consequence of the definition ( |2.13| ) for 0, we must have 



u ~ L to leading order. We therefore write the scaling relations 



and 



C„ (12) = ^F u (x) (4.15) 



C uu (12) = ^F uu (x). (4.16) 



The equations of motion developed in |T7| can be generalized to the post-gaussian case and 



produce the following relations between F u and F u 
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F u (0) 
+FJx 



7<> 



l-P)+2g 



WEZ 
p^/T—p 

2gf 



1 + v / W ? (l + v^T 1 



T\2 



(4.17) 



and 



4 



too + ^k(o) 

7T 



(4.18) 



For (7 = (|4.17|) simplifies to the gaussian form found in [17]. The quantity F u (0) enters 
into these equations and an analysis of ( |4.17| ) at x = gives 



^(0) 



(4.19) 



We then find that 



uJ, 



% 



1 - - 

7T 



>o, 



(4.20) 



which is a necessary condition for stability and is expected from the definitions ( f4.12|) and 
( |4.16| ). Equations ( [4.171 ) an d ( 4.18|) explicitly show how correlations in the u field are slaved 
to those of the order parameter. We have demonstrated here that the theory is consistent 
and that the fluctuations remain of (D(L~ 2 ). 



V. NUMERICAL ANALYSIS OF THE NON-LINEAR EIGENVALUE PROBLEM 



A. The post-gaussian theory without fluctuations 



The coupled non-linear equations ( pM6p , ( |3.27| ) and ( |3.33| ) compose an eigenvalue prob- 

(2) 

lem for the eigenvalues \i and A 2 which must be solved numerically. The eigenvalues are 
selected by matching the small-x behaviour, given by ( |3.36 ) and ( |3.37| ), onto the behaviour 
at large x, equations (|2.34j) and (|3.35|) . A fourth-order Runge-Kutta integrator is used to 
integrate ( |3.27| ) and (|3.33|) from near the origin (x = 0.001) into the large-x, asymptotic 
regime. Matching onto the proper large- x behaviour for T is the prime factor determining 



28 



//; the value of A% controls how well condition ( |3.35| ) is satisfied. The techniques used here 
are very similar to those used in previous studies 0. 

We have examined the 0(2) model without fluctuations in two and three spatial di- 
mensions. Table | contains the results for the eigenvalues fi and . The upper bound 
A% < 1/2 is satisfied both here, and later when we include fluctuations. The order- 
parameter auto-correlation exponent A can be computed using ( |2.32| ) once fi is known. The 
values for A obtained from the post-gaussian theory are presented in Table n, along with 



results from the gaussian theory ||17|| . The gaussian theory is in excellent agreement with 
simulations || of the 0(2) model in two dimensions, which give the value A = 1.171. The 
post-gaussian theory is in worse agreement with this simulation result. We will see below 
that the inclusion of fluctuations improves matters slightly. 

While the post-gaussian theory decreases A significantly, the form of the order-parameter 
scaling function T changes only slightly from the gaussian theory. The function T is plot- 
ted in Fig. [I] for the gaussian and post-gaussian theories in two dimensions. The minor 
difference between the two theories is reassuring because the gaussian theory is already in 
good agreement with simulations on this point [IJ. The functional form of / has the same 
qualitative features as T. The quantity g is shown in Fig. |2|. A key observation is that the 
first post-gaussian correction measured by g(x) is small for all x. 

The function C^p. (|3.18| ) measuring correlations in the square of the order-parameter 
field can be calculated from / and g. The physical, positive divergence in C^2 at small-x, 
which is seen in the gaussian theory, is now rendered negative. This occurs in two and 
three dimensions. One does not have this problem in the post-gaussian scalar theory j3|. 

(2) 

Superficially, this unphysical result is a consequence of A 2 not satisfying the lower bound 
( |3.45| ). More careful consideration indicates that the root of the problem lies in the method 
we chose to select P[m]. Our truncation of the expansion of P[m] in Hermite functional 
polynomials ( |3.1|) ignores terms in P[fh] that are 0(g 2 ). To remedy this we use a corrected 
form for C^z, 
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C^ = F 2 (l + g^J -1, (5.1) 

that differs from the previous form fl3.18| ) only at 0(g 2 ) and is manifestly positive near the 
origin. The BBS plot of 1/ (C^,a + 1) against 1 — T using the corrected form ( |5.1| ) is shown in 
Fig. |3]for the 0(2) model in two and three dimensions. The post-gaussian theory is in better 



agreement with the simulation results [III] than the gaussian theory. The two-dimensional 
post-gaussian results are in very good agreement with simulation data for 1 — T > 0.4, 
which are intermediate to large distances, x > 1. Unlike the gaussian theory, the shape 
of the BBS plot in the post-gaussian theory depends on spatial dimensionality. On this 
point the post-gaussian theory exhibits the same trends that are seen in the simulations - as 
the dimensionality increases from two to three the BBS plot approaches the gaussian result. 
These observations are strong evidence that post-gaussian statistics provide a more accurate 
description of the statistics governing the auxiliary field fh than do gaussian statistics. 

Fig. |] shows the scaling form g(x) ( |3.50| ) for the defect-defect correlations for both the 
post-gaussian and the gaussian theories in two dimensions. Simulation results |D| are also 
shown. While there still is a divergence in g(x) at small-x in the post-gaussian theory it is 
weaker than, and of opposite sign to, the divergence occuring in the gaussian theory. The 
relative weakness of the divergence is a consequence of the small value K 2 = 0.0018 in the 
post-gaussian theory, compared to K 2 = —0.5 in the gaussian theory. We see that the use of 
post-gaussian statistics does much to eliminate the the unphysical divergence in g(x), even 
in the absence of fluctuations. 



B. The post-gaussian theory including fluctuations 

The purpose of including fluctuations is four-fold. First, we want to render the auxiliary 
field correlation function analytic for small x. Second, we would like to completely eliminate 
the unphysical divergence at small-x in g(x). Third, we want to improve the agreement 
between the value of A obtained from the post-gaussian theory and the value seen in sim- 
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(2) 

ulations. Fourth, we would like to see if the inclusion of fluctuations increases A 2 ' so that 
( P-45 ) is satisfied. We choose B to have the form ( 2.13 ). In this eigenvalue problem ujq is 
selected using Q4.7|) , which guarantees that K 2 = L 2 = 0. We then solve the non-linear 
eigenvalue problem posed by ( |3.16|) , ( }4.3|) and ( |4.5| ) using the same methods we used above 
for the unmodified post-gaussian theory. For d = 2 the solution of the eigenvalue problem 

(2) 

gives uj = —0.2511, while for d — 3 one has uj = 0.1290. The eigenvalues /i and A 2 are 
given in Table |. For d = 2 the eigenvalue A 2 is increased; however, it still does not satisfy 
the bound (|3.45|) . In three dimensions actually decreases slightly. Table || contains 
the values for A. The agreement between the value for A obtained from the post-gaussian 
theory for d = 2 and the value from simulations is slightly improved when fluctuations are 
added. Adding fluctuations for d = 3 actually decreases the value for A slightly, moving the 
post-gaussian result farther away from the gaussian result. We cannot comment on whether 
or not this represents an improvement since, to our knowledge, no simulation data for A for 
the three-dimensional 0(2) model exists. 

The scaling form T for order-parameter correlations is only slightly modified from the 
form shown in Fig. [I]. On this point, the post-gaussian theory seems less susceptible to the 
perturbations introduced by the fluctuations than the gaussian theory ||17|| . Since the bound 
( |3.45| ) is not satisfied we must again use ( |5.1| ) for C^p. when creating the BBS plot. The 
results are shown in Fig. |5j. It should be noted that fluctuations do not modify the shape 
of the BBS plot in the gaussian theory because the leading order dependence of C^i and T 
on / is unaffected by u, which is 0(L~ 2 ). Again, we see that the post-gaussian theory is in 
better qualitative agreement with the simulation data than the gaussian theory. However, it 
appears that the addition of fluctuations removes the dependency of the post-gaussian BBS 
plot on dimension. 

The result for g(x) in the fluctuation-modified theory is shown in Fig. |j. We see that the 
inclusion of fluctuations not only eliminates the divergence of g(x) at the origin, but that 
g(x) is in better agreement with the simulation results down to smaller values of x than it 
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was in the unmodified post-gaussian theory. When compared to the gaussian theory with 



fluctuations |I7| we see that the addition of post-gaussian corrections has the desired effect 
of reducing the magnitude of g(0) and thus producing better agreement with simulation 
results. 

VI. DISCUSSION 

The main achievement of this paper is the extension of the gaussian theory for the 0(2) 
model to allow for non-gaussian statistics. In this sense the theory is unique. It is not trivial 
that the scaling properties of the gaussian theory are retained. We have demonstrated that 
when the structural form of P[m] is modified to include a post-gaussian part we can still 
solve the (now double) eigenvalue problem and find the fixed-point solution for P[m]. Thus 
we see a structure emerging, related to the nature of the fixed-point, that encourages us 
to believe that we can systematically change the structural form for P[fh] via (|3.ip and 
obtain improved results for the scaled quantities. Although the new function g is small, the 
approach we are attempting here is non-perturbative. 

In all areas examined, save the determination of the exponent A, the first post-gaussian 
approximation is in better agreement with the simulation results than the gaussian theory. 
Scaling results that are already in good agreement with simulations, such as T and the 
large scaled-distance behaviour of g(x) are essentially unmodified by the addition of post- 
gaussian corrections. On the other hand in the BBS plot (Fig. |3|), where there exists a 



large discrepancy between simulations [[LI]] and the gaussian theory, the addition of post- 
gaussian terms qualitatively improves matters. Unlike the gaussian theory, the BBS plot in 
the post-gaussian theory without fluctuations has a dependence on spatial dimensionality 
which shows the same trends seen in the simulations. In addition, the first post-gaussian 
approximation nearly removes the divergence in g(x) at small x. 

The development here suggests some avenues for improving the post-gaussian treatment. 
For now we have used the ad hoc formula which is correct to O(g), to obtain meaningful 
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results for the BBS plot. In order to properly treat these corrections at 0(g 2 ) we should go to 
the next order in the post-gaussian approximation sequence. In the current approximation 
the value for A is decreased from the gaussian result, and one hopes that the systematic 
inclusion of terms at 0(g 2 ) will raise A and provide better agreement with simulations. It 
is apparent that within this approximation scheme the value for A is converging slower than 
the results for the scaling functions, and may even be experiencing some type of oscillatory 
behaviour. 

It is also a non-trivial matter that we were able to incorporate fluctuations into the post- 
gaussian theory in a consistent manner. We have been able to successfully render both / and 
g more analytic and thus eliminate the divergence in g(x). The addition of post-gaussian 
terms to the fluctuation-modified gaussian theory has brought the value for g(0) into better 
agreement with simulations and thus generally improved the agreement at small-x. The 
fluctuations also slightly increase the value for A. We would also like to go beyond the 
the simplest post-gaussian fluctuation theory, represented by (|2.13| ). To do this it would 
seem that one needs to consider u as a post-gaussian field, coupled to the post-gaussian 
m Another interesting effect of adding fluctuations to the post-gaussian theory is the 
elimination of the dependence of the BBS plot on spatial dimensionality. Whether this is by 
accident, or is due to some deeper structure in the theory is presently unknown. It would 
be interesting to see if this effect remains if one includes 0(g 2 ) terms in P[m]. 

An obvious next step would be to examine the 0(3) model using post-gaussian statistics. 
The 0(3) model is interesting to study because of it's application to ferromagnetic materials 
and the role it plays in the description of monopole defects in nematic liquid crystals and 
cosmology |p5|| . There are some interesting new aspects associated with higher n |f2~4} . 



In examining the 0(2) model we have shown how the post-gaussian theory generalizes to 
a situation with continuous symmetry and solves some of the problems specific to the n = 2 
case. In developing the post-gaussian theory as a systematic and calculable extension of the 
gaussian theory we have also established the role that the gaussian theory plays as a zeroth- 
order approximation to the true statistics. There are many aspects of this process that 
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suggest a deeper underlying structure in the theory, including the interesting interactions 
between the post-gaussian corrections and the fluctuations. 
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APPENDIX: 

In this appendix we outline the method used to evaluate the post-gaussian average oc- 
curing in the defect-defect correlation function ^(r, t) ( |3.49| ). We have 

g dd (r,t) = (5[m(l))S[m(2)) de^V^l)] det[V 2 m(2)]> 

= e il .., n e i ,.. J „(A 1 (l) ■ • • D in (l)D h (2) ■ ■ -D jn (2)), (Al) 



with a sum over repeated indices and 



,0771,(1) 



^(1) = «K(1)]^ (A2) 



Using the definitions (|3.5| - |3.8|) and ( |3~9D the post-gaussian average in (|A1|) is 



(D h (l) ■ • • A„(l)%(2) • • -D jn (2)) = (D h (l) ■ ..D in {l)D h {2) ■ ■ ■ D jn (2)) 



+ Y: [ dld2A 2 (12)( [D il (l)---D in (l)D jl (2)..-D jn (2))) . (A3) 

k% J \6m k (l)6m k (2) f Q 

The first term on the right-hand side of (|A3|) is the original gaussian average, and was 



computed in JTj]. We focus on the term involving A 2 , which factors into the following 



product of n averages over the gaussian distributed scalar mj's: 
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From I|14| we know that 



(A 9 (i)i^(2)> 



h 



5j A + 



dh h 
dr r 



r~ r~ 

<*q Jq 



(A5) 



with 



2ir dr ' 



(A6) 



In (|A4j) the factor involving the integral is evaluated by operating with the functional deriva- 
tives, integrating by parts, applying the resulting delta functions and calculating the remain- 
ing gaussian integrals. We find 



dld2A 2 (12) 



5 2 



\5m k (l)5m k (2) 



[A fc (l)^ fc (2)] 



h 



dh h 
dr r 



with 



(A7) 



7r \ or or 



(A8) 



Thus the complete expression for g d d is 



g dd (r,t) = n\[-\ — + A 
\r I or 



(A9) 



with the post-gaussian terms contained in 



- (h 

k=l \ ' 



dh _h 

dr r 



l k Jk 



n 



h 



dh h 
dr r 



h q r 3q \ (A10) 



Due to symmetry, the only non-zero contributions to A come from terms containing either 
zero or exactly two factors of f. Mindful of this, we evaluate A and obtain 

ra-2 



A = n\ 



h 



dh h , ^hdh 
— -+ n-1 -— 

dr r r dr 



(All) 



The post-gaussian formula for g^d is then 

_ | (hV' 1 (dh dh\ (~Y~ 2 - — 

dd ' n \ r ) 1 dr dr j U U \r ) r dr' 



(A12) 
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TABLES 

(2) 

TABLE I. Values for the eigenvuales \i (top section) and A 2 (bottom section) for the 
post-gaussian theory, both with and without fluctuations. 





uj = 




d = 2 


0.8380 


0.7597 


d = 3 


0.5326 


0.5525 


d = 2 


0.05861 


0.08092 


d = 3 


0.06364 


0.05392 



TABLE II. Auto-correlation exponent A from the gaussian theory |l7j (top section) and from 
the post-gaussian theory (bottom section), both with and without fluctuations. 





LO = 


wo/0 


d = 2 


1.172 


1.269 


d = 3 


1.618 


1.655 


d = 2 


1.063 


1.092 


d = 3 


1.525 


1.513 
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FIGURES 




12 3 4 



X 

FIG. 1. Scaling form F(x) for the order-parameter correlation function in two dimensions. At 
x = 2 the upper curve is the post-gaussian theory without fluctuations and the lower curve is the 
gaussian theory without fluctuations. 
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FIG. 2. The auxiliary function g(x) for the post-gaussian theory without fluctuations in two 
dimensions. 
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FIG. 3. BBS plot for the gaussian and post-gaussian theories without fluctuations. At 
1 — T = 0.2 the upper solid curve is the post-gaussian result for d = 2, the middle curve is 
for d = 3 and the lower curve is the gaussian result. The solid circles are the simulation data for 



d = 2 and the open circles are the simulation data for d = 3 [11]. 
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FIG. 4. Scaling function g{x) for the defect-defect correlations in two dimensions. At x = 0.1, 



from bottom to top, the solid curves represent: the gaussian theory with ujq ^ [17]; the 
post-gaussian theory with loq ^ 0; the post-gaussian theory without fluctuations (diverging neg- 
atively); the gaussian theory without fluctuations (diverging positively) |Q. The dots represent 
the simulation data flTsl. 
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FIG. 5. BBS plot for the post-gaussian theory including fluctuations. At 1 — T = 0.5 the lower 
solid curve is the gaussian theory and the upper two curves, which are virtually indistinguishable, 
are the post-gaussian theory with loq ^ for d = 2 and d = 3. The solid circles are the simulation 



data for d = 2 and the open circles are the simulation data for d = 3 [11| 
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